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Abstract 

In this paper, we have applied the Meshless method based compactly supported radial basis 
function collocation for obtaining the numerical solution of unsteady gas equation. The unsteady 
gas equation is a second order non-linear two-point boundary value ordinary differential equation 
on the semi-infinite domain, with a boundary condition in the infinite. The compactly supported 
radial basis function collocation method reduces the solution of the equation to the solution of a 
system of algebraic equation, also, we compare the results of this work with some results. It is 
found that our results agree well with those by the numerical method, which verifies the validity 
of the present work. 

Keywords: Unsteady gas equation. Compact support radial basis functions. Unsteady flow. 
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1. Introduction 


1.1. Unsteady gas problem 

In the study of the unsteady flow of gas through a semi infinite porous medium [1] initially 
filled with gas at a uniform pressure Pq > 0, at time f = 0, the pressure at the outflow face is 
suddenly reduced from Pq to Pj >0 (Pi = 0 is the case of diffusion into a vacuum) and is, 
thereafter, maintained at this lower pressure. The unsteady isothermal flow of gas is described 
by a non-linear partial differential equation 

, , dP 

vV)-2A —, (1) 

where the constant A is given by the properties of the medium.In the one dimensional medium 
extending from z = 0 to z = oo, this reduces to 


d_ 

dz 


dP 



( 2 ) 
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with the boundary conditions 


P(z,0) - Po,0 < z < O 0 -, 
P(0,t) = Pi(<Po),0<t<oo. 


(3) 

(4) 


To obtain a similarity solution, Authors [2] introduced the new independent variable 


Vf V4Po’ 

and the dimension-free dependent variable y, defined by 

y(x) = -(1 - — 

a Pi 


0 


(5) 


(6) 


P-, 


where a = 1 - In terms of the new variable, the problem takes the form (unsteady gas 
equation) 


d^y 2x dy 

—^ +-^ = 0 , 

dx^ yjl - ay(x) dx 

x>0, O^a^l. 


(7) 

( 8 ) 


The typical boundary conditions imposed by the physical properties are 

y(0)=l, y(cx,) = 0. (9) 

A substantial amount of numerical and analytical work has invested so far [1, 3] in this model. 
The main reason of this interest is that the approximation can be used for many engineering pur¬ 
pose. As stated before, the problem of Eq. (8) was handled by Kidder[l] where a perturbation 
technique is carried out to include terms of the second order. Wazwaz [4] solved this equation 
non-linearly by modifying the decomposition method and Fade approximation. Parand [7-9] 
also applied the Lagrangian method, generalized Laguerre polynomials and Bessel collocation 
method for solving unsteady gas equation. Rad and Parand [11] proposed an Analytical solu¬ 
tion Unsteady Gas equation by Homotopy perturbation method (HPM). Kazem [14] applied the 
radial basis function (RBF) collocation method and Aslam noor [15] applied the variational iter¬ 
ation method (VIM) for solving non-linear this equation, khan [16] propose a new approach to 
solve unsteady gas equation, he applied the modified Laplace decomposition method (MLDM) 
coupled with Fade approximation to compute a series solution of unsteady flow of gas through 
a porous medium. Upadhyay and Rai [17] using the Legendre wavelet collocation method (yL- 
WCM) to solve this equation. 


1.2. CSRBF 

Many problems in science and engineering arise in infinite and semi-infinite domains. Dif¬ 
ferent numerical methods have been proposed for solving problems on various domains such as 
FEM[35, 36], FDM[31, 35] and Spectral[8, 13] methods and meshfree method[5, 12]. The use of 
the RBF is an one of the popular meshfree method for solving the differential equations [32, 34]. 
For many years the global radial basis functions such as Gaussian, Multi quadric. Thin plate 
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spline. Inverse multiqudric and etc was used [6, 10, 33]. These functions are globally supported 
and generate a system of equations with ill-condition full matrix.To convert the ill-condition ma¬ 
trix to a well-condition matrix, CSRBFs can be used instead of global RBFs[22]. CSRBFs can 
convert the global scheme into a local one with banded matrices. Which makes the RBF method 
more feasible for solving large-scale problem [23]. 


1.2.1. Wendland’sfunctions 

The most popular family of CSRBF are Wendland functions. This function introduced by 
Holger Wendland in 1995 [24, 25]. Wendland starts with the truncated power function 0;(r) = 
(1 - which be strictly positive definite and radial on for / > L|J H- 1 , and then he walks 
through dimension by repeatedly applying the operator I. 

Definition [26] with 0;(r) = (1 - r)^ we define 

<f>s,k - ( 10 ) 


it turns out that the functions f^.k are all supported on [0,1]. 

Theorem 1 [26] The function <ps,k are strictly positive definite (SPD) and radial on M' and are 
of the form 


fsAr) 


Ps,kir) re [0,1], 
0 r > 1, 


( 11 ) 


with a univariate polynomial of degree L|J + 3k + 1. Moreover, 6 C^^{R) are unique 
up to a constant factor, and the polynomial degree is minimal for given space dimension s and 
smoothness 2k [26]. Wendland gave recursive formulas for the functions f^.k for all s,k. We 
instead list the explicit formulas of [27]. 

Theorem 2 [26] The function f^.k, k - 0, 1, 2, 3, have form 


fsfi = (1 - r)[, (12) 

=(l-r)';‘[(/+l)r+l], (13) 

= (1 - r)';2[(/2 +41 + 3)r^ + (31 + 6)r + 3], (14) 

4>s,3 = (1 - -I- + 231 + 15)r^ -H (6f + 361 + 45^ (15) 

-H(15/H-45)r-H 15], (16) 


where I -'i^\+k+I, and the symbol = denotes equality up to a multiplicative positive constant. 
The case k - 0 follows directly from the definition, application of the definition for the case 
k - \ yields 


fs,\(r) 



(m(r) 


t(l - ty^dt — 


^00 

Jr 

r 


tfMdt 
f(i - tydt 


1 

(i+m + 2) 


(l-r)'^‘[(/+l)r+l]. 


(17) 

(18) 
(19) 


where the compact support of fi reduces the improper integral to a definite integral which can 
be evaluated using integration by parts. The other two cases are obtained similarly by repeated 
application of I. [26] We showed the most of wendland functions in Table 1 . 
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Table 1: Wendland’s compactly supported radial function for various choices of k and s=3. 

<ps,k smoothness SPD 

03.oW = (l-r)2 C" m3 

03.iW = (l-Ot(4r+l) C2 m3 

<Px 2 ir) = (1 - r)® (35r3 + 18r + 3) M3 

033(r) = (l-r)^(32r3+25r3 + 8r+l) C® M3 

03’4(r) = (l-r)|°(429r4+45Or3+21Or3 + 5Or + 5) C* M3 

<Px 5 {r) = (1 - r)f (2048r3 + 2691 + 1644r3 + 566r^ + 108r + 9) C'° M3 


1.2.2. Wu’sfunctions 

In [28] we find another way to construct strictly positive definite radial function compact 
support. Wu starts with the function 

0(r) = (l-r3);,/eiV, (20) 

which in itself is not positive definite, however, Wu then uses convolution to construct another 
function that is strictly positive definite and radial on M 

Hr) = Jjl - hlH - (2r - (21) 

then, he constructed functions by ’’dimension walk” using the D operator. 

Definition 2 With <pi(r) = ((1 - * (1 - .3)^)(2r) we define 

(pk.i - D'‘cf>i. (22) 

The function are strictly positive definite and radial on M® for i < 2A: + 1. are polynomials of 
degree 41 - 2A: + 1 on their support and in in the interior of the support. On the boundary 

the smoothness increases to [26] 
a number of Wu’s function showed in Table 2. 

Table 2: Wu’s compactly supported radial function for various choice of k and 1=3. 

fs.k smoothness SPD 

0(),3W = (l-Ol(5 + 35r+lOlr3 + 147r3 + lOlr4 + 35r3+5r® C® M3 

<Pl 3 (r)^(l-r)l(6 + 36r + S2r^+ 12P+30r*+5r^) H M3 

i*2!3W4(l-r)3(8 + 40r + 48r3 + 35r3+5r4) M3 

i*3!3W = (l-r)t:(16 + 29r + 20r3 + 5r3) C® M3 


The Oscillator [26, 29] and Buhman [30] functions are the other kind of CSRBFs can be 
showed in Tables 3 and 4. 
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Table 3: Oscillator’s compactly supported radial function for various choice of k. 


(pk 

smoothness 

SPD 

(pi{r) = (1 - r)\{l +Ar- 15r^) 

C 2 


02 (r) 4(1- r)6(3 -H 18r H- 3r^ - \92r^) 



03 (r) 4(1- r)*(15 -H 120r -H 210r2 - 840r^ - 3465r4) 

c® 



Table 4: Buhman’s compactly supported radial function. 

smoothness SPD 


<p\{r) - 12 r^logr 
^ 2 {r) — 2 r‘* log r - 

03 (r) = r* - —r® + 


■2lr^ + 32r^' 


74 16 3 

—r H- r 

2 3 

1024 


r-2 - 378r4 + 


99 

04(r) = 


132r® ■ 


9216 


- 12r2 + 1 

2 r 2 +i 

. 1024 7 84 , , 

——r2 ——r + 1 

11264^, , 396 


-r 2 


35 


-r2 


198r^ 


35 


+ 1 


C 2 




2. CSRBF method 

2.1. Interpolation by CSRBFs 

The one-dimensional function y{x) to be interpolated or approximated can be represented by 
an CSRBF as 

N 

y(x) y„(x) = ^itpiix) = (23) 

(-1 

where 


(f>i(x) = (pi — 

- 2 Cill^ 
r,,. 

(24) 

(t>^(x) = [4>\ix), (plix), • • • , (Pn(x)], 

(25) 


■■ ,^Nf, 

(26) 



(27) 


X is the input, is the local support domain and ^,s are the set of coefficients to be determined. 
By using the local support domain, we mapped the domain of problem to CSRBF local domain. 
By choosing N interpolate points {xj, j - 1,2, • ■ • , A^) in domain: 

N 

yj^Y = 1,2, • • • , A^). (28) 

i=l 

To summarize the discussion on the coefficients matrix, we define 

AE = Y, (29) 
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where : 


= [yi,y 2 , 


(30) 

= [(l)^(vi),<l>^(V 2 ),-- 

• ,<l>^(x^)]' 

(31) 

'(^i(vi) <^2(21:1) ■ ■ 

■ Mm)' 


<^l(-^ 2 ) (plixi) ■ ■ 

■ Mm) 

(32) 

,(Pi{Xn) (plix^) ■ ■ 

(pNiM- 



Note that (pi{xj) - (p( 
achieved. 


Il-^i - XjW 

r<j 


), by solving the system (29), the unknown coefficients will be 


2.2. Solving the transfonned Model 

In this section, by defining u{x) as a below form and transforming ^ and ^ in terms of it. 


u(x) = 1 - ay(x) 

(33) 

dy du 
dx dx 

(34) 

d^y d^u _j 

dx^ dx^ 

(35) 


the problem of unsteady gas in a semi-infinite porous medium (8) takes the following form : 


d^u 2xa du 
dx^ yfu dx 


v>0,0<a'< 1, 


and the conditions change to 

m(0) = 1-0", m(oo) = 1. 

Now we approximate ^ and ^ as 


du du„ 
dx dx 

d^u d^Un 
dx^ dx^ 



By using integral operation u{x) is obtained as 


(36) 

(37) 


(38) 

(39) 



(40) 

/^X 

u„(x) - u(0) = V 1 (f>iir)dr, 

(41) 

^X 

(x) = (V 1 (l>i{r)dr) -H (1 - a), 
tli -'0 

(42) 
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By substituting (38), (39), and (40) in (36), we define residual function 


Res(x) - 



2x du„ 
■\fu^ dx 


(43) 


Now, by using N - I interpolate points plus a condition (37) we can solve the set of 

equations and consequently, the coefficients will be obtained; 


^^h(^cxd) 1 ■ 

. Collocation point are chosen by : 

;•= l,2,...,iV-l, (45) 

where p is a arbitrary parameter. Here, we choose Loo = 5 which satisfies y(Loo) < e with e as 
a small positive value. The parameter g and local support domain that must be selected by the 
user. But here, by the meaning of residual function, we try to minimize ||/?ei(x)|p by choosing a 
good local support domain and g. We define ||/?ei(x)|p as 

nL^ ^ T T 

||/?ei(x)||^ = I Res^{x)dx djRes^i-^Sj + (46) 

Jo ^22 


where 






, ; = 0, l,...,m. 


T’m+i(sj = 0, ; = 0, l,...,m. 


(47) 

(48) 


Pm+\{.s) is the {m + l)th-order Legendre polynomial. The Figure 1. show the minimum of 
||f?ei(x)|p which is obtained with for cases ofWendland^ s, VTM 33 , BuhmanA, and Oscillatorj,. 
The result of this section can be summarized in the following algorithm for the BVP; 

L(x,y(x), ^, ^) = 0, y(0)^a, y(oo) ^ b. (49) 

dx dx^ 

Algorithm The algorithm works in the following manner: 

(1) Choose N center points from domain [0, 

(2) Approximate ^ as the from ^ 

(3) Obtain y{x) by using defined integral operation I^(h{x)) — h(t)dt in the form u^ix) - 
Ijti X"" 4>iWt + a 

t i2 

(4) Substitiute un(x), and into the main problem and create residual function Res(x). 

(5) Substitiute collocation points {Xj}y~Q into the Res(x), along with a boundary condition y(La,) = 
b and create N equations. 

( 6 ) Solve the N equations with N unknown coefficients {^,)^[ and find the numerical solution. 
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Figure 1: The graph of residual error versus local support domain when N = 30, a = 0.5 and xi =1.5 

3. Results and discussion 

In this section, we compare the applied result of CSRBF with RBF methods [14] and yLWCM 
[17]. 

In the physical observation of the unsteady gas problem, y'(0) has an important issue [1]. Table 5 
and 6. presents a comparison between the values of y'(0) and y{x) for a - 0.5 and obtained by 
the CSRBF and, RBFs and yLWCM. 

In addition, the value of y'(0) obtained by CSRBF methods for n - 20,30 and various values 
of a are reported in Table 7. 

The stability of the CSRBF scheme depends on the local support domain r^. An important 
unsolved problem is to find a approach to determine the optimal size of Also, The condition 
number grows with N and g for fixed values of local support domain. Figure 2. and 3. shows the 
condition number of matrix A versus local support domain and 

Figure 2. shows that by decreasing the r^, an increase is seen in the condition number of 
matrix A and the method become more unstable. In general, the smaller the value of the 
higher percentage of zero entries in matrix A, so condition number matrix A is a little and with 
increases the r^j, decrease percentage of zero entries in matrix A and matrix A will be unstable. A 
percentage of Zeros and condition number of CSRBFs in Matrix A 30 X 30 based on different values 
of local support domain illustrated in Table 8. 

To find the critical values of , we use the residual error which can see in Figure 1. Figure 

4. displays the y{x) for a - 0.5 in Unsteady gas equation. 

in Figure 5. displays the convergence rates of the CSRBFs method for Eq. (8) for some value 
of N for a - 0.5. These figures illustrate the convergence rate of the method. 
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Table 5: Comparison of >'(;r) for a = 0.5 with g = 1.5, = 30. 


X 

Wendlands^s = 2.5 

W 113.3 = 2.3 

Oscillator^ = 3.4 

Buhman^ = 2.8 

RBF.G 

RBF.T 

yLWCM 

0.1 

0.88136427 

0.88136428 

0.88136409 

0.88136468 

0.88139802 

0.88137298 

0.88147552 

0.2 

0.76582809 

0.76582823 

0.76582774 

0.76582852 

0.76588029 

0.76579834 

0.76661101 

0.3 

0.65599963 

0.65999935 

0.65599915 

0.65600036 

0.65606928 

0.65590792 

0.65727018 

0.4 

0.55389758 

0.55389797 

0.55389693 

0.55389814 

0.55399431 

0.55375746 

0.55567752 

0.5 

0.46094112 

0.46094164 

0.46094037 

0.46094191 

0.46107202 

0.46078383 

0.46325882 

0.6 

0.37797968 

0.37798027 

0.37797895 

0.37798037 

0.37814380 

0.37783286 

0.38078076 

0.7 

0.30535020 

0.30535087 

0.30534921 

0.30535091 

0.30554174 

0.30522806 

0.30853930 

0.8 

0.24295205 

0.24295279 

0.24295099 

0.24295286 

0.24316554 

0.24285666 

0.24644872 

0.9 

0.19033171 

0.19033246 

0.19033052 

0.19033235 

0.19056419 

0.19025895 

0.19409149 

1.0 

0.14677064 

0.14677143 

0.14677694 

0.14677139 

0.14702048 

0.14671538 

0.15877555 

1.2 

0.08313306 

0.08313387 

0.08313169 

0.08313370 

0.08341060 

0.08310227 

- 

1.4 

0.04404326 

0.04404412 

0.04404186 

0.04404399 

0.04433767 

0.04402800 

- 

1.8 

0.01002990 

0.01003085 

0.01002839 

0.01003049 

0.01034206 

0.01002778 

- 

2.2 

0.00170828 

0.00170909 

0.00170673 

0.00170889 

0.00202559 

0.00170941 

- 

2.6 

0.00021351 

0.00021440 

0.00021195 

0.00021413 

0.00053177 

0.00021524 

- 

3.0 

0.00001696 

0.00001786 

0.00001540 

0.00001758 

0.00033532 

0.00001878 

- 


1.3702e-08 

5.843 le-9 

2.4820e-08 

1.4394e-08 

8.52e-07 

- 

- 


Table 6: Comparison of initial slope >’'(0) for a = 0.5 with g = \.5, N = 30. 


X Wendlands^s r^j = 2.5 

Wu3,3 = 2.3 

Oscillator^ r^^ = 3.4 

Buhmariix r^j = 2.8 

RBF.G 

RBF.T 

yLWCM 

y’(0) -1.191796 

-1.191806 

-1.191800 

-1.191768 

-1.191498 

-1.191243 

-1.199258 



Buhman . Wendland - Wu-- Oscillator 


Figure 2: The graph of condition number versus local support domain when N = 30, a = 0.5 and ^ = 1.5 
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Buhman . Wendland - Wu-- Oscillator 


Figure 3: The graph of condition number versus ^ parameter when N = 30, a = 0.5 and = 3 



X 


Figure 4: solution of >'(x) when N = 30. 
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Table 7: CSRBFs solutions of initial slope >’'(0) with N = 20,30 


a 

Wendlands^s 

WU2,3 

Oscillator-^ 

Bithman^ 

N-20 

N-30 

N-20 

N-30 

N-20 

N-30 

N-20 

N-30 

0.25 

-1.15838196 

-1.15658845 

-1.15868018 

-1.15652004 

-1.16347841 

-1.15658199 

-1.15914956 

-1.15651901 

0.50 

-1.19674806 

-1.19179615 

-1.19388602 

-1.19180634 

-1.19468231 

-1.19180040 

-1.19476755 

-1.19176821 

0.75 

-1.23833094 

-1.23998794 

-1.23307537 

-1.24033353 

-1.23372756 

-1.23980995 

-1.23493307 

-1.23984689 


Table 8: A percentage of Zero and condition number in Matrix Asoxso- 



Zero 

Wendland^^ 

WUt,2 

Biihman^ 

Oscillator^ 

1.0 

64.44% 

7.789055e04 

4.532194e04 

6.667960e03 

6.196000e03 

1.5 

50.45% 

5.005580e06 

4.490489e05 

4.762594e04 

1.137094e05 

2.0 

38.55% 

1.070687e08 

2.277201e06 

2.017929e05 

9.349735e05 

2.5 

28.00% 

1.284523e09 

8.029373e06 

6.255295e05 

4.865744e06 

3.0 

19.33% 

5.416534e09 

2.248322e07 

1.581147e06 

1.88648le07 

3.5 

12 .00% 

6.652775el0 

5.373891e07 

3.469944e06 

5.974023e07 

4.0 

06.33% 

4.377016el0 

1.142168e08 

6.860756e06 

1.642163e08 

4.5 

02 .11% 

8.799029el0 

2.230297e08 

1.252852e07 

3.99027 le08 

5.0 

00 .00% 

1.338415el2 

4.055797e08 

2.148696e07 

8.539759e08 



Figure 5: The graph of residual error versus N when = 3, a = 0.5 and xi =1.5 


4. Conclusion 

The fundamental goal of this paper was the construct an approximation to the solution of 
unsteady gas equation, a set of CSRBFs with these properties were proposed for providing an 
effective but simple way to improve the convergence rate. A comparison was made among the 
solutions of [14, 17] and this work, the absolute error ||Rei||^ were obtained. This paper has 
provided an acceptable method for Unsteady Gas equation. It was also confirmed by logarithmic 
figures of residual function that this method has an exponential convergence rate. Additionally, 
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high convergence rate and good accuracy are obtained by the proposed method using relatively 
low numbers of collocate points. 
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